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GOPDARP SPACE FLIGHT CENTER 
Greenbelt, Maryland 


ELECTROSTATIC HEAT fejX INSTABILITIES 


INTRODUCTION 

The purpose of this report is to investigate plasma physics of the Electrostatic Ion Cyclotron (EIC) 
and Ion Acoustic (IA) modes in the presence of a heat flux. The initial analysis will be carried out 
mdre generally to include a current in the plasma as well as a heat flux. This will facilitate com- 
parisons with the results of Kindel and Kennel (1) who investigated these modes driven by a current 

/,'/ 1 O 

I ■. v-> v > 


We wish to find a dispersion relation for oo(k) parameterized by the ratio of electron to ion temper- 
ature and also the growth rate, T, as a function of the wavenumber k. We then wish to specialize 

■ v * ti 

these results for the special case of marginal stability for which T(k) = 0 and find the threshold 

drift speed above which the plasma is unstable to the EIC and IA modes. Finally we will find ex- 

ff % _ 

pressions for these quantities for the limits — 1, and ^ ^ . 

\\ ^ 

A self-consistent distribution function will be derived in the next section using Chapman-Enskog 
theory to represent an electron distribution including a heat flux characterized by a thermal 
conduction speed, v 0 and also an electric current characterized by a drift speed, u. The ions will 
be represented by an isotropic Maxwellian distribution. In the third section this distribution will 

«a 

be used to find the dispersion equation and growth rate for the low frequency^ CtM< JL*,* n i c ; 
electrostatic modes, In section IV marginal stability will be assumed while the contribution of 
the electron current will be dropped and the dispersion equation simplified for the two modes, 

EIC or IA. Here an expression will also be given for the minimum conduction speed above which 
the plasma becomes unstable. The limits ^ 1 and ^ 1 

will be considered and simplified expressions for the dispersion equation, growth rates, and mini- 
mum critical conduction speed derived, The final section will then consist of a short discussion 

j ' 

of the results. Oi 
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THE ELECTRON DISTRIBUTION FUNCTION 

In Chapman-Enskog theory It is assumed that the particle collision frequency, v, Is very large so 
that f (the distribution function) can be expended in a series; 

f * Ao * y, A, * . . . 

m +* f * *■ f . + . . , 


where the A,’s are all assumed to be of the same order so; 

It. I *lfj *>U a l • • • 

In addition is is assumed that all derivatives of the A ( ’s are of the same order as A ( 
derivation we follow Tanenbauin. (2) From ilk equations (A.5.20) and (A, 5. 21): 


itself, In this 


where: 


r * = - -£_r 

T I -V <V I. 




neglecting any pressure gradient and C a AT - u. where U. is the peculiar velocity and 


fa a Af T>c • ^ so off-f 0 +f, in this approximation leads ^o 

the NavienStokes equation. Our distribution function is now; 


Tfc«A*/»L 






where f 0 is 


Now let us find some moments of this. In particular consider the heat flux; 
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We assume that there is a temperature gradient in the Z direction, parallel to the magnetic field, Then; 


where; 

Then: 


c*« 


since 


0 x i t 

V| «" * Art* % •*£ 


(9) 


C • /*T" Vk *nd 4x • ** % • 


( 10 ) 


so 
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P Of in C. n 
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and the distribution function in terms of the beat flux is: 


( 12 ) 


where: 


* s! ~Vm 


A » iHirA 'i 
» ■ 

Now let us calculate the first movement, the average velocity. We can see that fj is even in v x and 


v v while it is odd in v, x * v, 

y 

and <?>• <*;■> £ , 

It is therefore usefu 


since £ is directed only along the z-axis. Consequently 




then: 
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Then f can be re-written as: 


f«* a[i» a *.£*!!♦ C K€ " ** S) . 

Then, in general: ^ 

and integrating over : 

Consider the odd moments: 

<*r> * 3# fcw-*> «■****.. 

o 

* a 

where the first term in the integrand has vanished because of its symmetry. Using: 

e -*^ x . <£M^! 1 ^ 

where: 

OJ/n-O!! * 1*3 • 0 S * . , Xa/n-l) 

We see: 

<'<«»*’“’> * [ta«*0 1 ! - 3 Uk- 01 ll ** 

which yields the general result: 

(-K-0I! B • 

Then, the first moment is: 


04) 


(15) 


<*r>- -** 


w* it ■ 

<AQ„'> r O and since Al^ #l * _ aTV* 


Then: 


(16) 


(17) 


< aT„> ■ u. . 


(18) 

(19) 

( 20 ) 
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We can then interpret the parameters as follows; is the field aligned heat flux density 

(and is collision dominated) while u is the electron drift velocity associated with a net current 
■a ^ j 

flowing parallel to 0 .We note that the choices of u and are independent so setf/ng 0 

will reproduce the results of Kindel and Kennel while setting* * 0 will specialize to he'it flux in- 
stabilities which are the main concern here. 

To complete an analogy with the cuii'-jt instability it is useful to define one other quantity; This 
is the thermal conduction speed, . From equation (10): 


/ 


o , » - § . 

^ •• * ** m H ->•« •** * 


( 21 ) 


Replace by the mean free path \^and the temperature gradient yjjf by where L is the sys- 
tem’s scale size. <£> (| then becomes; 


9., * “ 3 ’l’ 1 [ * 


( 22 ) 


which is in the standard form for a heat flux: a speed multiplied by an energy flux density. The 

O' 3 

heat conduction speed can then be defined as; v 


*■« ■ 


(23) 


where the factor 3 is put in for later convenience. 

This definition can be justified on physical grounds as follows: Consider an ideal gas. The heat 
energy density is nc p T since heat conduction will assume constant pressure (mass motions caused 
by the pressure gradients are taken care of by the parameter u in the distribution function), The 
heat flux is then: 

g - -a r c (*tc f T) 


(^ c v T ) 


&' 


..... •*., »,*» Jj .'M , . 
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and for an ideal gas; <V *** while 

<P » - f aT c (i 


! '/c, 


so; 


9*-i/T c ^i a T W4b ° VC * 

Consequently the constant £ is not involved in /£ , Thus the parameter B can be written as: 

6 » < 24 > 
and this form will be useful later. 


We can further note that there is a limit to the magnitude of since, in the collision dominated 
regime assumed in the derivation, it can't exceed the product of the random electron flux (across a 
plane;perpendicular to the z-nxi.s) and the mean thermal electron speed: 


r 




( 25 ) 


which, in concert with equation (22) implies: 


i 


or: 


<X 4J 


( 26 ) 


and this is what is meant by flux limited. 


THE DISPERSION EQUATION 


For two charge species (electrons and ions) the dispersion equation is (3) : 

* O * 1 ♦ Xe ♦ : * ° 


( 27 ) 
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assuming: 


^ 9 % 4 


8 * 6* a 

*r « *n,t ■* aJV* 4* 


*^ A « at h 4 4 *r/ 


o a fevAfc 

JC -< /**!., C! 

“*<* -7C?~ 

a«< » «»* sr« H * “xir 


A 3 /vT • A^yt/*lA/r H . 

ions so it assumes, for 4 X + O 


This expression is solely for electrostatic perturbations so it assumes, for 4 

■, / 

a *r^cx/H4>. raersuftc . « 

M 11 MA6«tert« PA,essro*e >> -L - 

For the ions, the distribution function is taken to be an isotropic Maxwellian: 


where fflt is the ion thermal speed: 


fti* 


4 *« r < 

/mi 


Then, from Hasegawa* 3 *, equation (2.30): 


v»* * *& to T& 




where: 


- i {bpf 

- modified Bessel function 

= Plasma dispersion function (Fried & Conte* 4 *) 
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so; , 

(V 

For the electrons we will use the distribution discussed above: 

(** aJi* B^-C^Oa^-siO « 


i^ktl 4 4J|J ) 


( 31 ) 


( 32 ) 


where: 


a » * 

8 * ’ 




and: 




/w« 

* ” 

With these definitions: 


ft i. - u. 


aTI 


«x 




4'\ 

V3 


wartj 3 * 


( 33 ) 


where: 


X e ■» * ax,^ . 


Now for ion cyclotron and ion acoustic instabilities 1 so H 4 <«1 . Then the 

argument, X, of the J is appreciably different from O only for /*•*,»> X . However, 
there is little contribution to the integral from large /*e± due to the factor €" in f 9€ . Con- 
sequently, we will assume << 1 and with 3^, CX^di O (n 1) and 


for 1*0 and the electron term becomes .* 


«t4o 




o 

( 34 ) 


II'.' 


where; 






Now; 


jfe * - 4S ]] e t ' r * ?4 " i * 5 


(35) 


and the integrals over all have the form: 


i <■ 


H4I - 


* «*! . 


Performing the integrals over /Q . then gives; 




where: 


* 


^ t*0 = at.,, ^ -c,r € ? ♦* Jl . 


Now we define: 


"2.k(S) ’ TT- J 


K -V 

x e 

*-5 


<&* 
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(36) 


(37) 


(38) 


(39) 


so that 2.CS) is the plasma dispersion function of Fried and Conte (4 \ The electron’s con- 
tribution is then: 

fc. « A [2, tf.) ♦ &(* 2,1 f.> - <• 2, CU * I C40) 
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& * - • 

If 2 h is integrated by parts the following relations may be found; ^ 

e . * - i e.' 

u ■ -iS.K 

0 

2k- it 

and H/» -a (,1 ♦ 3 2/1 where * j[j* * 

Successive Iterations of the recursion relation gives; 

wiiile successive applications of the relation between 2 % and 2, gives 

2," s «Ur-l)2o + H y 

Combining these gives: 

* ja.i- i * $€. «•! 

i«) - resell) * J4, 

V?) -r 4 2. tj a *i »?«.♦* 

V*> + j o - ■ je. 



( 46 ) 


so X« can be re-written entirely in terms of 2, as: 

*€ • j $8&V [i 

' (C-/' - 

At this point is useftil to consider the magnitudes of B and f € in order to make approximations. 
We have already shown that! 

G< i. c^Tc 

since is flux Ur sited. Now consider t 

G - » Am, U. WV » U. 

where is the wave’s phase velocity parallel to the magnetic Held. For T */n x i 
we expect (following Kindel and Kennel) that m. a* r** * if u is taken as the drift speed for 
marginal stability. Similarly, in the same temperature regime we expect sff c ** 
so 

(t 

S « « 1 *>1 , 

When we actually have an expression for <*>(£)we must check to ensure that this limit is correct. 

In what follows, we will assume <*JL in order to Uiop higher orders of tothedisper- 
sion equation. We must keep in mind that this is strictly true only for r */r L >*> 1. and 
when we write the general result we will include these terms, We will, of course, also include them 
when Finding results in the t */t*c £ £ regime. Note that the dispersion equation looks 
increasingly like that of Kindel and Kennel as higher order terms in are dropped so there is 
some heuristic value in this approximation. The corollary to this is that our results should increas- 
ingly deviate from theirs as becomes large. ’ 


Now, keeping only terms which are explicitly of first order in & or ; 


( 47 ) 


The dispersion equation is: 

4 ? 0 

1 ♦ 1 « 4 - X i - O . 

(( 

Using equations (30) and (47), dividing by 2, this becomes: 

n 

i) 

o. 


(48) 


(49) 


Considering the ion contribution, T , the imaginary part of V , is alwsys positive and so pro- 
vides damping. ,The n * 0 contribution is Landau damping while for i ' , it Is referred to 

as ion-cyclotron damping. ,^J) 


Considering the electron contribution, ?0 so the sign of 

the sign of ^ 5 + which is the factor multiplying 2, *. 



\ 

*Tf - u. — wy 


C 


where 


depends on 





Consequently, if the drift velocity, u, is larger than the wave’s phase velocity then the negative 

rs 

electron dissipation can contribute to instability. In the same fashion, the heat conduction 
velocity, if sufficiently large, can lead to instability. The mode will become unstable if the veloc- 
ities, u and , are sufficiently large that the negative electron dissipation is larger than the 
damping from the ions. 

O o 

GENERAL CASE 

In general, if there are a number of ion species, X* must be replaced by T X i where the sum 
is over all the ion species. With this and the definition: 






Si ..an awf as 
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T m ^-C ‘1*00 


the dispersion equation becomes; 

IV. A C" T’ Tm(* 


+ jTc^» [l- 1 - C^sSr 5 ) Zj; 4 * ^X^"] 


whete; 


* a *Ki "T1 
A*i * 


li _ 1. f 

* *■ l 'Z 3 Ci J 

A * i* T; 


» x 

J' _ <• — 

" * 3 /*»• **9 AT** 


where the last term includes the higher order terms in ^ previously dropped. 


Now, assume the system is near stability so: 

6j * + j, y 


AaI * 




Then set u « in equation (50), separate into real and imaginary parts: 

D* &*^n a 40 ♦ ofc;U» Aj O 


and then, to order ZT ; 

A 

E>< - o 

T» b i 

1 ’ " ”HT" 

(Krall&Trivelpiece (S) , 1973). 


ufe&Mtr.. 


o 




Defining: 


2* « R«U«) 

2* - 3^ lO 


o 


gives the real dispersion equation (multiplying by A* X^): 

-aA-V. 

♦1 ♦ § a»C«aaf) 

- !(£)$£) [^ 4(W- 


(55) 


we can see that the conduction speed, a£ , introduces new terms in the dispersion equation that ’re 

Y * { i 

not seen in the case of current instabilities alone, These terms all contain the factor: 

(^n)m 

which vanishes if /Q vanishes but not when u vanishes, 

o 

From equation (50), 0 ; («*y , can be re-written as: 

where: C : - 

This gives the combination of drift speed and conduction speed for marginal stability, The terms 
in braces multiplying aQ are those dropped when $ € «JL 


(56) 


Using: 


* rrc 


5 * 


(5/7) 


'V) 
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the real imaginary parts of the dispersion equation for small growth rates can be re-written in 
terms of our normalized variables as: 

-IX Jt 


‘ r> 


!s.s« [ i - J.U (!.*- 4) w}_* O.-l) w.) 


(58) 


*♦*[&$ Wl)} • 3fcc 


(59) 


where; 


ft * $1 

*T-iSr 

lyn ® ^ 

« - «*■ 
Je’ /r« e 

arfm ££&£ 
% 

/» / \*4/ Te\fe 

f = ii. 

„ S3- 

r.oe * tirffi.e* 

ir* 

i» 




and the underlined terms can be neglected for | € « 1. . 


To find the marginal drift speed for instability equations (58) and (59) are solved simultaneously 
along with the requirement that u. be a minimum with respect to k^ and k n . Note that 

by equation (54) this is equivalent to requiring P =0 unless ^*^0**0 . That is, what we wish 
is to find, for a given set of physical conditions (T e , Tj, n e , etc.) the solution to the dispersion 
equation, eufk^), for which the combined drift speeds are minimum. 


Using equations (54) and (50) the growth rate is 

T » To* 

/i " “OoTTom 
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( 60 ) 




i 


where; 


Tof 




(60 


Bottom, IX % + 


* /* 


Z £ f (l-3j;) - 5j e (■'4-3)2«.'i -[i* is.\£ i>] <a\} 


(62) 


where the variables are defined identically to those used in equations (58) and (59) and equation 
(57) has been used. We also note that since Jf c X is a requirement the term in J c In equa- 
tion (62) can generally be dropped unless ^ is appreciably greater than unity, At this point it 
might be useful to re-iterate the assumptions implicit in equations (58) through (60): 

1) low frequency: <**n<*Jl.* 

2) Electrostatic modes only 

3) Chapman-Enskog approximation ** collision dominated plasma ** /'fi, </J^ e or < 1 

4) Heat flux parallel to § -field 

o 

5) Maxwellian ions 

6) iu*| . 

HEAT FLUX INSTABILITIES 

In what follows we will set u * 0 so there is no electron current and no net drift velocity, We will 
also assume only one ion species so n t * n„ and the sum over ion species will be removed. With 
■* 7**Te the equations now become: 

r - (63 

l! s.cv- 3 .*)-[i iq*, 13.3 j 
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( 64 ) 


3 *c rrtsTciT 






] 


As long as J C «X the major difference between these equations and the analogous equations 
for the current driven instability is that the conduction speed doesn’t occur in the arguments of the 
electron z-functions in the heat flux case while the drift speed does appear in these arguments in the 
current driven case* When ’ Vr i *>1 , «1 , and the equations for the two cases 

are identical. We can then directly take over the results of Kindel and Kennel, replacing the drift 
speed (their V ) by the heat conduction speed, At high temperature ratios then the mode is ion 


acoustic in nature with 


« 1 

■Si* 


. This is done below. 


On the other hand, when 1 , the critical conduction speed (the conduction speed 

above which the mode is unstable) becomes of order the electron thermal speed. 

Then the terms in Ji will be important and the nature of the instabilities will differ somewhat 
from the current driven case. As we shall see the instability takes on the nature of an ion cyclotron 
wave for X characterized by largely perpendicular propagation, "Vt.’O 1 ' 

This is true for both the current driven or heat flux driven instabilities but the details are some- 
what different. 


The major difference between the two cases occurs for 


. In the case of a heat 


flux driven instability the critical conduction speed becomes larger than the electron thermal speed. 

i) 

Since it is not possible to conduct heat faster than the thermal speed of the particles the plasma is 
stable with respect to this source of fi^e energy. To transport energy more rapidly than it is 

necessary to direct the electrons rather than let them diffuse: This is a current and gives rise to the 

e 

current driven instability Tor sufficiently large drift speeds. 

We will now summarize the results for t Vr 5 >•> j. (since the analysis is identical to that 
given in Kindel and Kennel), briefly consider how the results differ from the current driven case 

17 




when *t/n ,w 1 , and then present gmphs of the numerical results for both ranges. 

We will consider the ion acoustic approximation separately from the ion cyclotron case and show 
that the former has the lower critical drift speed in this temperature regime and so the mode is 
actually ion acoustic in nature, 


0 Ion Acoustic mode: Reaming to equations (63) and (64) we set k A * 0 so the wave is prop* 
agating solely along the/^eld lines. This implies thatu^ * 0 and since: 


I'm Co) 


AntinilAfio KnnAniA* 
lilt* v\|iiUMV?na> vvywiMvt 


m* o 
/nfo 


- Jf: * L1«-S^Z *(■?.} 


(65V 


& ' $A~ 


i» C: e 


r a 


1- | Je*(i-5«*) 



(66) 


where « "? * was used in the second equation, We expect a result which will roughly be 
cmJ *5 ^ (| C, which would then imply; 



o 
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2.CJ3* 


We will therefore assume $ e « \ ^ and use the approximations: 

-A* jW 

for the electron and ion plasma dispersion functions. The dispersion equation then becomes; 

a c » 


f 


(69) 


co « * r-Ttrfy— (70) 

j ^ , is the sound speed (k B is Boltzmann’s constant) which matches our ex 
pectations. The equation for the conduction speed is now: 


. Vi; 

where C. * 




To find the mininum conduction speed we differentiate this with respect to /d"y 


ui 

A,. 


equate the result to zero: 




(71) 

an d 

(72) 


which has the approximate solution for large Cj : 

■C * 


(73) 


where is the value of which minimizes X* Using this in equation (7 1) yields the 


minimum marginal conduction speed: 






(74) 


or for a hydrogen plasma; 


1 


Jw.s- + tJU T *V 
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(75) 





for the ion acoustic mode with r */n >> 1 , For large T «/n the critical drift speed is 

only a few times the ion thermal speed. 

In this temperature regime rf/Xi is generally large (from eq (70)) so ion Landau damping is 
weak. If the conduction speed exceeds the parallel phase velocity then the electron distribution 
will have a positive slope at a speed equal to so the electron Landau instability will cause the 
wave to grow. 


Now consider the growth rate in this limit. With the above approximations it becomes: 

.3* n a 

71 « * _ Ttt. g * 

(j?) “V3* Hf 1 


(76) 


it 


and since L *• j 5?. the last two terms in the denominator may be neglected, Then, using 

r. > * 


equation (70): 




r - 3s] ti- a ] 


(a) 


(81) 




where x* Cl + V 


1-. 


® Electrostatic Ion Cyclotron Mode; Following Kindel and Kennel we note that for ion cyclo* 
tron waves *“1 in this temperature limit so we can neglect w. *•«, . Further we can 

use 3* - m 3 C y "O X. for “11 the ion terms since ^ (l is small for cyclotron waves. For 
temperature ratios near unity 40^-Jl.i but as T Vn increases to will increase to keep damp- 
ing from the n = 1 mode small. If to were to get close 2S2, then the n * +2 damping will become 
important. Thus we will have oO x £ JL; . Then Kindel and Kennel find the minimum 
critical conduction speed to be; 

^ l(L V 

| 82) 
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or* for a hydrogen plasma; 


£7 * 


( 83 ) 


The perpendicular wavenumber Is roughly given by; 

(%£*■) ~ 1 - »'*:(%) 


(84) 


and: 


9 fu «r £ 1 


(85) 


As T */r t increases the wave changes to more nearly parallel propagation. We also see this in 
comparing equations (83) and (75); the ion acoustic wave has a smaller critical drift velocity then 
the ion cyclotron. The mode becomes increasingly more like an ion acoustic wave as the tempera- 
ture ratio increases, Because of this it is not useful to calculate a growth rate in the ion cyclotron 
limit. 


T« 

T* 


1 : 


From the numerical solution of equations (63) and (64) we see that when <* 1 then 

3 c **•5' while 3« ' w *1 . Consequently, even though they don’t make a large contribu- 


tion, the terms in in equation (63) should not be dropped. We will however drop the terms 


m m 

in 3 in the denominator of equation (64). 


0 \> 


Then, keeping only the largest terms, we have: 




( 86 ) 


^ ’■ ahJkjimAxsi^SSta^S^bt£stt&. 
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Jwll* 

#■* 


Cli- 





(87) 


Now, from the equation for the growth rate we see that we must have ^ fitf) 

in order that the electron terms contribute to growth rather than damping. Consequently we will 
expand for small' J € and drop terms of order compared to f c . Then 

equation (86) is: 


In the ion acoustic case in the high temperature limit this changes equation (70) to; 


( 88 ) 


/W.\* g/ 

> X* * <3 J. S. 


p 



and since S c > 0 this means that the phase velocity would be less than expected when this 

term is ignored. 


© Ion Acoustic mode; 

As before we set k^* 0 solv^* 0 to get: 


o 


*«■ f. I . 


(90) 

(91) 


We note that the corresponding limit in the current driven case has equation (90) without the last 

, X\ a 

term. There is no solution to that equation for positive \ since the maximum value of 

2.* is about .57 at Sc )• S . In fact there is no solution for ^ O 

for r «/Ti t 3.r. 
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The inclusion of the extra term in equation (90) wlU allow solutions for lower values of r *//t than 
this. Equation (91) is identical to (71) and so leads to the same minimization of eq (72). Solving 
(72) numerically for ^ * 1 yields; 

fim * *72 
- .017 


Substituting these values into equation (90) gives: 


W 


•su *■ 


-.62, 


i! 


This violates the initial assumption that k is real and w complex. This same situation also occurs for 
smaller values of t Vt<. . 

The identical problem appears in the case of a current driven instability and is discussed by Kindel 
and Kennel. The resolution is that the smallest critical value of is zero (so the real frequency is 
also putatively zero). However, a small increase in the drift velocity above the critical value results 
in a shift to w * Wfi (the ion plasma frequency) and to .In practice 

only ion oscillations near . will be observed. We shall also see that near the mode is 

basically Ion Cyclotron rather than Ion Acoustic which is why the ion acoustic approximate gives 

\?r"' 

a poor result, 

We see that as ^ decreases the critical conduction speed increases from a few times the ion thermal 
speed to nearly the electron thermal speed. When t Vtc -* 1 the wave’s phase velocity is 
roughly the same as the ion thermal speed for conduction speeds not greatly exceeding the critical 
value. For example, when n -2 the most unstable ion acoustic wave has *** f • 
so ion Landau damping is strong and a large electron drift is required for instability. If "^in- 
creases to large values then •»■> and ion damping can be neglected. This leads to the 
Buneman (6 * instability, 


o 
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Since the critical value of £ lies near unity no simplifications can be easily made in the equation 
for the growth rate. 


<B> Electrostatic Ion Cyclotron Wave; 
For this wave ^•Aj.-sOL so in general! 






a- Ml; 


V> X, 


(92) 


will be assumed and then equation (69) can be used to approximate the ion dispersion function. 
Further we wiU see that *-i so we can assume 

\o« ^ 1 

and drop this term from the equation. Also we will again expand 2* in the small limit. 

Thus equation (88) becomes*. 

A 1C _ V" - 'Jlf ? 

I* 4 ' U - L_ C4-/nJt.i 3^*1* . C93) 

Now we require that while still satisfying equation (92), This again requires that k /( , 

be small. Then using the addition rules of Bessel functions equation (93) can be re-written as: 


!♦% -acM ■■ 


Jk&a 

U -Ai 



while the equation for the critical conduction speed is: 



(94) 


(95) 


where ^1*0 ■* T* • * (96) 

“Solving” equation (94) for to (ignoring the fact that and are functions of c o would give 

<a) ** -*.: U+ fe) (97) 
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where A anti we see that the Inclusion of the term 4jf tends to decrease L 

putting o) closer to n,. It is not particularly worth while to carry the analysis further for this 
case. The numerical solution of equations (63) and (64) (with the requirement that><£ be a mini- 
mum) is given on the following pages and contrasted with the solution to the case of a current 
driven instability and as demonstrated above there Is no significant difference until 1 , 


As can be seem from these results and those of Kindel and Kennel, the minimum critical conduction 
speed is smaller for the ion cyclotron mode than for the ion acoustic. Consequently the instability 
takes on the characteristics of an ion cyclotron wave as the temperature decreases to unity. We alsc^ 
note that as decreases to unity increases so there is some temperature ratio below which no 
physical heat flux can drive the plasma unstable; the threshold to too large. 

Let us now find an expression for the growth rate for the ion cyclotron mode using the same set of 
approximations. Froth equations (60-62): 




(98) 


or; 


V 






-’WJ. 

1 f* 


(99) 


Now, (Mi)< 1. for all A. since T7V « 1 and £ * TU * Thus ^e ter ws in the sum 

t? 

are smaller than unity except perhaps for n = 1 and since for large n the sum behaves as 

? 

T tS" and so the contribution to the sum is small for large n. An order of magnitude ap- 
proxlmation to the sum is thus given approximately by the n = 1 term; 


where & is defined by equation (97) and is small, Since the solution here is not very different from 
the current driven case we can drop the last term in equation (94) to get an approximate solution 

for 6 : c. 6 

* ~ 1- & * T;/r«. 


where a is defined by equation (96), The approximation to the sum is then: 

~(i-a*V)*- 

and since \<i ( »*) \ < X. while ** X. for this case the sum is within an order of 

magnitude of unity, Since it is multiplied by Jr and 1 we see that the last two 

terms in the denominator of equation (99) can be neglected since they are less than unity. Con* 
sequently, the growth rate can be re-written as: 


JL 

co 


9 P . T* (o IL 

<T" '"TV.UO 


5“ ** ■*« tA; 


(100) 


for ^ ^ and h given by equation (97). 

Tc * ^ 

On the following pages are given graphs of the dispersion equation and growth rate as well as the 
minimum critical drift speed, and M; all as functions of % 
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Figure 1. Minimum Critical Drift Velocity 
as a Function of Te/Ti 











DISCUSSION 

The major point to note is that, unlike the current driven instability, the heat flux instability will 
not “turn on” when^r < 1 since the requirement on the minimum critical heat flux is that it be 
unphysically large (ie,,f cm > 1), The only way in which electrons can transfer energ^Tastcr than 
the electron thermal speed is by directed (rather than random) motion. This is a current and gives 
the current instabilities. 


The second point to note is that the graphs given here for the heat flux instability do not markedly 

differ from those Kindel and Kennel have produced for the current instability, This has already 

Te 

been mentioned in the analysis for^r » 1, 

From the graphs it can be seen that the results here differ from Kindel and Kennel’s by 10-20% 

Tt» t« _ 

nearer * 1 and arc essentially identical when is greater than 3 to 4. This means that the two 

Te 

instabilities can be “combined” when ^ la greater than unity. That is, the fact that u and v c 
don’t enter into the real and imaginary parts of the dispersion equation in an identical manner is 

*T' 

Irrelevant for^T > 1. Thus we can take the critical drift velocity curves and replace either u (in 

Te 

Kindel and Kennel’s curves) or v c (in our curves) by the combination u + v c for^r > I. For 

Te 

example, the minimum critical velocity for the instability to occur with - 3,0 is about ,12. 

Thus we can say that if the sum of the heat flux drift velocity and the current drift velocity: 

J)(( v c + u > . 1 2 

then the instability will occur, Thus, the existence of a heat flux can cause a current carrying plas- 
ma to be unstable to the EIC or IA modes even if it were not unstable with the current alone (and 
vice versa), In a solar flare, for example, the conditions which lead to a current also will give rise 
to a heat flux and the combination make the plasma less stable than either alone, 


ThOhirdjJoint of interest is that the two modes, EIC and IA, exist as separate entities with the 

\\ yip a Te 

loweisminl'mum critical velocity for the EIC mode when^r < 8,5 while the IA mode has the lower 


minimum critical velocity for > 8,5. The EIC mode is characterized by > l and the IA 

k Tg y.;. 

mode by ^ < 1Q“ S (from the numerical solutions). Consequently when is somewhat less 
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than or greater than 8,5 the two modes will both occur, one propagating parallel to the field and 
the other perpendicular to it, if the actual drift velocity is somewhat greater than the minimum 
critical value, For example, if the drift velocity is roughly 30% above the minimum critical value 


in solar flares depends on the presence of the EIC mode 7 this gives a wider possible rangeiin tem- 
perature over which it can occur if the drift velocity is only slightly above the minimum value, 


Some results not found in the graphs or the analysis but which are seen in the numerical results 

are that some parameters can be changed by large amounts with correspondingly small changes in 

v c . In the EIC mode calculations it was found that the minimum critical drift velocity is only 

Te T e 

weakly dependant on the value of tjj-, For example, at * 2.0, a change of 20% in 

only a 2% change in v c and this is a typical result. Thus the growth rate of the EIC mode is not 

strongly peaked about one value of ; when v c slightly exceeds v cm we would expect that the 

k 

destabilized wave will contain a spread in Thus the wave will propagate in a cone around 
the magnetic field and the angular thickness of the cone will be appreciable compared to the angle 
the side of the cone makes with the field line. 

For the IA mode, the preferred direction of propagation is very nearly exactly parallel to the 
magnetic field with k A /k t t < 10“*, On the other hand, although v cm depends strongly onjj, 
if we write / 1 as: 

then we find that changes in/ C yc and as large as a factor of 2 (but which leave/ { constant to 

six significant figures) change v c only in the sixth significant figure. We expect then that the IA 
mode, when unstable, will propagate with a wide range in w. 

One further point can be made. For Te/Ti » 1 we found that the IA mode was unstable with 
lowest v c and dispersion equation given by equation (70), If the term k* * 2 m tlle denomina- 
tor is ignored this further simplifies to; 


k 

produced 

■I I 


<** s * kf | C $ J 


on 


,, M /Te 
Jim V 

We see that the computer generated graph of S i m under these conditions duplicates this result to 
within about 10%, If we use this in equation (71) we find that, for^sr-» 1: 


'em 


v te 

and this depends on temperature only very weakly. 


=“.0174 


SUMMARY 

We have investigated the plasma physics of the Electrostatic Ion cyclotron and Ion Acoustic modes 
as driven by a heat flux, We find that no physical heat flux can make the plasma become unstable 

fr. ° nn 

to these modes when < 1 . For ^ in the range 1 to 4 our results differ from the case of a current 

Te 

driven instability by less than 20% while for^i" greater than this the two situations are virtually 
identical. We can consequently use the curve for the minimum critical drift speed with this speed 
replaced by the sum (v c + u) m where v c is the electron conduction speed and u is the current drift 
speed. 


T'f* 

The general dispersion equations have been given as well as simplified versions for^~ I and 
Te 

jj-» 1. The exact equations have been solved numerically and graphs of the results given. 

If 

In ar. appendix the first order correction term in is given for the Ion Acoustic mode. 

k n 
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Appendix A,* First order correction term in£*- for the Ion Acoustic Mode, 

*i 1 


f- 


The lack of information on comes from the approximation p, » 0, tet us consider the first 
order correction in p,; 

r m (M|) -e^« I m (iUj) 


and: 


e** t *l-^ 
I 


so to the first order in p. : 

Then, with ? c f e « I equation (63) becomes: 

!+2kf, x’a, - ff j WZj, Iff) + (Z R Cf,) - !4Z R <r, - f cyc )] ) + (f c - f t ) z R (f t ) 
and the same approximations in equation (64) yield: 

f c «f e { l+qe-fi 3 + te J (J-p,(l-J4eMc(2^-l))]| 

Now, if we assume fj and fj - £ cyc » I and £ e « 1 (so we can neglect the last term in equation 
(A5)): 

1 +2k ll *de *(jrT)-jjr 
using equation (69). This can be re-written as: 

i 
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(A2) 

(A3) 


(A4) 

(A5) 

(A6) 


\»i »cyc / 


(A7) 


o> 

rsr 2 > 


i^uxi - l 4u^\ (^) 


(A8) 


or: 


['-^y 


so we can clearly see the correction term in 

k H 

“Solving” this for <*) gives: 


(A8) 
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a‘ 


kiiCsl 


1 + 2k| | Xj, 


: / Ai — 2 


«fe) ' ft) ' (a 

\ ft. 



(A9) 


co 


so we can see (since rr>> 1) that the correction term reduces the frequency below the standard Ion 




ft 


Acoustic result of equation (70). 

In terms of the variables used in the computer solution equation (A8) can be written as; 


or, solving for co; 


to* 


Te 

Or 

1 

|j ~‘A w 2 v 2 

(A10) 

2 (v- 

iVi 

2 x u 2 v 3 \ 

(All) 


2 A 

Te/Ti / 


where: fj * uv, v « — and x * 1 + 2kJ| . We can estimate the last term from equation (73): f j m 2 st 




w 2 *■ 


/ v - 1 \ / 2 x In 2 CA 

\v-2/Y Te/Ti ) 


(A 12) 


and for large > x ** 2. Using this value: 


■# 

Te 

4 In 2Ci 

Ti 

Te/Ti 

10 

3.2 

20 

1.8 

30 

1.3 

40 

50 

1.0 

,83 

100 

,45 


Tg 

so there is no solution for « until > 40. This gives a rough estimate of where the approxima- 
tions are valid, In this range v » 1 so roughly: 

co=*£“E (A 13) 

Where k is the factor in the last set of parentheses in (Al 2). 

Te 

As an example, for * 50: 
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